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> '. Abstract 

I The two-photon correlation of the Ught pulse emitted from a sonolumines- 

I cence bubble is discussed. It is shown that several important information 

, about the mechanism of light emission, such as the time-scale and the shape 
P-t'. 

of the emission region could be obtained from the HBT interferometry. We 

c , 

, also argue that such a measurement may serve to reject one of the two cur- 



rently suggested emission mechanisms, i.e., thermal process versus dynamical 
Casimir effect. 

The sonoluminescence process converts the acoustic energy in a fluid medium into a short 
light pulse emitted from the interior of a collapsing small cavitating gas bubble. Since the 
discovery of the technique for trapping a single cavitating bubble by standing acoustic wave 
many remarkable properties have been revealed The spectrum of emitted light 

is very wide, extending from the visible to the ultraviolet regions. The conversion process 
requires some mechanism of extraordinary concentration of the energy and is apparently 
related to keeping the sphericity of the collapsing bubble while its radius shrinks more than 
one order of magnitude. An important fact is that the light emission takes place within 



a very short period of time compared to the typical scale of hydrodynamic motion of the 
bubble. The sensitivity of the radiation power and spectrum to physical parameters such as 
temperature, pressure, the amplitude of acoustic drive and the gas composition is the most 
intriguing aspect of the process and no apparent reason nor any well-defined systematics are 
known. In particular, the emission mechanism of light is still controversial. Some authors 
attribute the light emission to the quantum-electrodynamic vacuum property based on the 
dynamical Casimir effect 0-^. Others consider that the thermal process |1T0|-[12| like black- 



body radiation should be the natural explanation for the process. Non-equilibrium atomic 
collision process could also be a strong candidate [l^. In any case, the gas dynamics inside 



the bubble |]14|,|T5|, in particular the shock wave formation, seems to play an essential role 
T6| . However, experimental information on the dynamics of the gas inside the bubble is not 



available at present. Due to the short time scale and to the smallness of the emission region, 
the precise measurement of these geometric and dynamical information is quite difficult. 
Usually, the time dependence of the bubble radius is measured using the Mie scattering 
process of laser beams 0,|T3]- Nevertheless, since we have no information on the properties 
of the gas during the implosion phase, it is not obvious that the scattered laser amplitude 
is really measuring the bubble surface. On the other hand, besides the knowledge of time 
elapsed in the process, it would also be desirable to determine the shape of the emission 
region. 



The two-photon interferometry, initially proposed by Hanbury-Brown and Twiss 



for measuring stellar sizes, later found its application in the analysis of high-energy particle 
production processes and is nowadays widely employed in relativistic heavy-ion collisions 



|, |21|| . The basic principle is associated to the Bose- Einstein statistics obeyed by the 
identical particles involved in the process and to the chaoticity of the emission mechanism. If 
the emitting source has no additional dynamical correlation, then the two-particle correlation 
is directly related to the geometrical size of the source. Several authors |2^J2^ have discussed 
the effect of the source dynamics on pion interferometry. In this note, we apply this well- 
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established method to obtain some important information on the dynamics of the gas inside 
the bubble while it emits the light pulse. It is interesting to note that the application of 
HBT interferometry to the sonoluminescence bubble {R ~ 10~^m) lays just in between the 
stellar [R ~ 10^°m) and the high-energy physics [R ~ 10~^^m) scales. 

The HBT interferometry method consists in measuring two light quanta in coincidence. 
Let ^2(^1 5 ^2) be the probability for simultaneously detecting two photons with wave- vectors 
ki and ^2 and Piiki) the single-photon probability. The correlation function C{ki, is then 
defined as 

c{kM)^ ^^fi'^^^ • (1) 

If the emitting source is chaotic and static, the correlation function is simply related to the 
Fourier transform of its space-time distribution, so that from the measurement of C(A;i, ^2) it 
is possible to determine the source geometry. In general situations, though, the relationship 
between the correlation function and the source geometry is more complex, being not possible 
to determine uniquely its geometry only from the knowledge of C{ki, k2)- Nevertheless, once 
we have a good guess about the shape and the time development of the emitting source, we 
could suggest an appropriate parametrization for the geometry of the process, and the HBT 
measurements can be used to determine these parameters. 

It should be emphasized that the chaoticity of the emitting source, which manifests as 
random phases of emitted signal, plays a crucial role in relating the correlation function 
to the geometry of the source. In contrast to this, if the emission process is coherent (for 



example, a laser source), it is well-known p0| , p4| , p5| that no HBT correlation would be 
observed, that is 

c{k^M) = 1. 

As an immediate consequence, we can distinguish these two extreme scenarios by a precise 
HBT measurement. As mentioned before, the mechanisms of light emission proposed so- 
far can be classified into two categories: One based on rather conventional atomic process 



of thermal origin, and the other due to the dynamical Casimir effect. In the latter case, a 
coherent burst of light would be produced, whereas for the former scenario we expect chaotic 
emissions. In both cases the single-photon spectra are similar to that of the black-body 
radiation, but an HBT analysis would allow for clearly differentiating these two emission 
mechanisms, shedding some light on the most intriguing feature of the sonoluminescence 
phenomenon. 

Let us now investigate some examples of non-trivial HBT correlation functions (C 7^ 1). 
For this, we assume that light quanta emitted from different space-time points has no extra 
correlation besides that due to the Bose-Einstein statistics. This is the case when, for exam- 
ple, the emission process has its origin in atomic collisional or bremsstrahlung mechanisms. 
We further assume that the gas inside the bubble be locally in thermal equilibrium. This 
does not necessarily mean that the emitted photons are also in thermal equilibrium with the 
gas. For simplicity, we here only treat spherically symmetric sources 0. 

Since the time scale of the hydrodynamic motion is considered to be few orders of mag- 
nitude greater than the time scale of the emitting process (~ lOps), we first consider that 
the source size remains constant during the light emission. Thus, the time and space de- 
pendences is factorized. In Table 1, we show the analytic expressions of the correlation 
function for some typical parametrization of the source density p{r,t). The first three cases 
(A, B and C) refer, respectively, to a Gaussian type of source, to a spherical shell and to a 
homogeneous sphere, all of them having a Gaussian lapse of time. The fourth example (D) 
corresponds to an exponential spatial distribution shining constantly within an interval of 
time r. 

It could well be possible that the emission region is not static, in which case the time 
dependence of the source size should be considered. One may, for instance, imagine that 
the radiating source is the spherical domain behind the expanding shock front formed at 
the center []TT],|TD[. In any case, the fluid velocity is certainly much smaller than the speed 



of hght so that any dynamical effect due to the fluid motion could be completely neglected. 
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An emitting source with these characteristics is represented by the case E. Throughout this 
study we have taken h = 1. 

TABLE 1 





Form of the Source 


C(ki,k2) - 1 


Case A 






Case B 


6{r - R) e-*'/2r2 


g-(A..)V2 [sin(gi?)/(gi?)]2/2 


Case C 


Q{R-r) e-*'/2r2 


gg-(A..)V2| [cos(gi?) - sm{qR)/{qR)]/{qR)Y/2 


Case D 


e-r/R 0(3^2 _ ^2) 


'sm{AiuV3T)/AujVSTf {l + q^R^y^/2 


Case E 


0(i?t-r)e-*'/"'e(t) 


9|/|V(8/u'), 

/ = -t^ [(1 + ijz+)W{z+) - (1 - iJZ-)W{z-)] - 2/i 



Here, q = ki — k2 , Au = lji—lj2 = c{ki — k2), fi = Rt q, z"^ = (Aa;±i?g)r/2 and W{z) = 
e~^^erfc {—iz). R is the spatial extension parameter and r is the time span parameter of the 
source and R is the velocity of the shock wave. The case A was studied by Trentalange and 
Pandey p6|, following the prescription given in However, as pointed out by Slotta and 
Heinz [E8|, the current conservation should be correctly taken into account, which results in 



a slight deviation from eq.(6) of Ref. [^]. Note that in the first four cases (A-D) above, the 
correlation function are written in the form 



C{k,,k2) = l + -T{Au;)^q) 



(2) 



which is a consequence of factorized form of the source. Therefore, to determine the time 
span parameter r, it is convenient to plot the data for fixed q. If the source is really 
factorized in space and time, then plots of log(C — 1) vs. (Acj)^ for different pairs of {ki, ^2) 
should generate a set of identical curves just shifted from each other for different values of 
q = ki — k2 ■ In Fig.l, we show examples of such plots, with r = Ips, R = Ifim in the 
blue-light domain (fci ~ ^2 — 4 x 10^ m~^). The solid lines refer to the case A, the filled 
circles to the case D and dotted lines to the case E. As can be seen in this figure, if the 
time span is Gaussian and factorized from the space dependence, then these curves should 



be parallel straight lines whose slopes would give the value of the parameter r. If the time 
span is not Gaussian, then the curves are not straight lines, but similar to those for the 
case D. However, even in this case, the first derivatives of the corresponding curves with 
respect to (Acj)^ at the origin provide an estimate of the parameter r. The information 
about the source spatial distribution, ^{q), is determined by means of the intersection of 
these lines with the abscissa. For the case E of Fig.l, the velocity of expansion was taken 
to be i? = 2 X 10~^c, where c is the speed of light. If the velocity is smaller than this value, 
the lines corresponding to different g's come closer to one another. This happens because, 
for smaller velocities, the effective emission region becomes smaller and, consequently, the 
correlation function becomes broader and slow-varying with q. Note that for the case E, the 
lines are not parallel, reflecting the non-factorized emission source. 

In practice, such an analysis might be limited by experimental conditions. In particular, 
to get a meaningful result for the r parameter, Au should be measured within a resolution of 
the order of 1/ps, which corresponds to measuring the photon energy itself within precision 
of 6uj/uj < 10~^. This might not be easily achieved. However, even if the energy resolution is 
not high enough, the spatial factor in eq. can independently be analyzed if the source is 
factorized. Being so, the HBT correlation function of two photons of approximately the same 
energy, under a poor energy resolution (within 6uj), would lead to an angular correlation 
given by 

Since the multiplicative factor {T^Su)) is independent of q, it can easily be eliminated from 
C near the origin, q = 0. The necessity of the above renormalization of the correlation 
function at the origin as a consequence of a poor energy resolution of HBT measurements 
has been pointed out in The function in (^) behaves quadratically in q around 
g = 0. Namely, 
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1 - Ry + ■■■ 

1 - RY/^ + 
1 - Ry/5 + ■ ■ 



Gaussian (A), 
Shell (B), 
Sphere (C), 



1 - 4i?2g2 H , Exponential (D), 

so that we can determine the parameter it! by the curvature of ^{q) at the origin. By defining 



K, 



q=0 



we have 



^ K,/2, Gaussian (A), 
^3k/2, Shell (B), 

R = ( 

5k/2, Sphere (C), 
Exponential (D). 

To distinguish among the shapes of the correlation functions refiecting different source den- 
sity distributions, we would need to know the behavior of ^{q) in a wider range in q. To 
stress the differences among the four cases above, let us introduce the variable X defined as 



X=^q. 



Then, by definition, $ behaves as 



$ - 1 - + 



near X — Q. In Fig.2, we compare the behavior of $ as a function of X for the four 
cases. The continuous line corresponds to the Gaussian spatial distribution, the triangles 
to the sphere, the squares to the spherical shell and, finally, the circles correspond to the 
exponential density distribution. As can be seen in this figure, the differences among these 
curves are not striking near the origin, but if the data are precise enough (~ 3 order of 
magnitudes) in a sufficiently wide range of X, we may determine the shape of the source 
function p. 



As already mentioned, if the emitting process were coherent, the behavior of the corre- 
lation function would be entirely different from the ones just discussed. As anticipated in 
the beginning, a precise measurement of two-photon correlation function of the light quanta 
emitted by a sonoluminescence bubble would allow for distinguishing between chaotic and 
coherent emission mechanisms. However, we should stress that, even in the case of a chaotic 
source, a poor energy-resolution experiment could lead to a result similar to the one expected 
in the case of a coherent source, due to the factor exp (— r^Ao;^). Therefore, a very high 
energy-resolution experiment is required to clearly differentiate between these two opposite 
scenarios. 

This work has been partially supported by FAPESP (Contract 95/4635-0), CNPq and 
FAFERJ. One of us, (T.K.) expresses his gratitude to J. Rafelski, I.Scott and H-T. Elze for 
stimulating discussion on sonoluminescence phenomena and encouragements. Discussions 
with C.E. Aguiar, L. Pimentel and P.A.Nussenzweig are gratefully acknowledged. 



8 



REFERENCES 

[1] D.F. Gaitan and L.A. Crum, J. Acoust. Soc. Am., Supple. 1 87, S141 (1990), D.F. 
Gaitan, L.A. Crum, C.C. Curch and R.A.Roy, J. Acoust. Soc. Am. 91, 3166 (1992). 

[2] Sonoluminescence, L.A. Crum, Physics Today, Sep. (1994). 

[3] B.R Barber and S. Putterman, Phys Rev. Lett. 69, 3839 (1992). 

[4] B.P. Barber, C.C. Wu, R. Lofstedt, P.H. Roberts and S.Putterman, Phys. Rev. Lett. 
72, 1380 (1994). 

[5] R. Heller, K.Weninger, S. Putterman and B.P. Barber, Science 266 248 (1994). 

[6] K. Weninger, S. Putterman and B.P. Barber, Phys. Rev. E54, R2205 (1996). 

[7] Casimir Energy for a Spehrical Cavity in a Dielectric: Toward a Model for Sonolumi- 
nescence?, K.A. Milton, |hep-th/951009Tl . 

[8] C. Eberlein, Phys. Rev. A53, 2772 (1996). 

[9] Sonoluminescence and the Heimlich Effect, A. Chodos, |iep-ph / 9604368 . 
[10] R.Lofstedt, B.P. Barber and S.J. Putterman, Phys. Fluids A5, 2911 (1993). 
[11] W.C. Moss, D.B. Clarke, J.W. White and D.A. Young, Phys. Lett. A211, 69 (1996). 
[12] Sonoluminescence: Sound into Light, Setli Putterman, Scientific AmericanFev. (1995). 
[13] L.Frommhold and A.A. Atchley, Phys, Rev. Lett. 73 2883 (1994) 
[14] I. Scott, H-Th. Elze, T. Kodama and J. Rafelski, submitted to Phys. Rev. Lett. 
[15] M.C. Chu, Phys. Rev. Lett. 76, 4632 (1996). 
[16] C.C. Wu and P.H. Roberts, Phys. Rev. Lett. 70,3424 (1993). 
[17] W.J. Lentz, A.A. Atchley and D.F. Gaitan, Applied Optics 34, 2648 (1995) 
[18] R. Hanbury Brown and R.Q. Twiss, Phil. Mag. Ser. 7, Vol. 45, 663 (1954). 

9 



[19] G. Goldhaber, S. Goldhaber, W. Lee and A. Pais, Phys. Rev. 120, 300 (1960). 

[20] See for example, M.Gyulassy, S.K. Kauffmann and L.W.Wilson, Phys. Rev. D20, 2267 
(1979). 

[21] Proceedings of the Int. Workshop on Correlations and Multiparticle Production 
(CAMP/LESIP IV), ed. by M. Pluemer, S. Rahaand R. M. Weiner, World Scientific 
(1991). 

[22] Y. Hama and S. Padula, Phys. Rev. D37,3237 (1988). 

[23] S. Pratt, Phys. Rev. D33, 1314 (1986), K.Kolehmainen and M. Gyulassy, Phys. Lett. 
B180, 203 (1986) 

[24] R.J. Glauber, Quantum Theory of Coherence in Quantum Optics, edited by S.M. Kay 
and A. Maitland, Academic Press (1970). 

[25] C-Y. Wong, Introduction to High-Energy Heavy-IonCoUisions, p431. World Scinetific, 
1994. 

[26] S. Trentalange and S.U. Pandey, J. Acoust. Soc. Am. 99, 2439 (1996). 
[27] D. Neuhauser, Phys. Lett. B182, 289 (1986). 

[28] Photon Spin and the Shape of the Two-Photon Correlation Function, C. Slotta and U. 
Heinz, |nucl-pli/96100T6| . 



10 



FIGURE CAPTIONS 

FIG.l Correlation functions plotted as functions of Ao; for varios values of q. The solid 
lines refer to the case A, filled circles are for the case D and dotted lines are for the 
case E. The first two cases are examples of factorized sources, whereas the lines for 
the case E are not parallel. 

FIG. 2 Geometrical form factor ^{q) plotted as functions oi X — ^ —^^-^^Q- The sohd 
curve corresponds to the Gaussian source, the triangles to the spherical source, the 
squares to the spherical shell and the circles to the exponential distribution. 
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